1 



cn 
o 
o 

(N 



> 

O 
O 

a^ 

O 
m 
o 

Oh: 



X 



^^/(iV) Gauge Theories Near Tc* 

B. Lucini'^t^ m. Teper^, U. Wenger^* 
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We study the deconfinement phase transition in SU{N) gauge theories for 7V=2,3,4,6,8. The transition is first 
order for A'^ > 3, with the strength increasing as A*' increases. We extrapolate Td ^/a to the continuum limit 
for each A'^, and observe a rapid approach to the large A'^ limit. As A'^ increases the phase transition becomes 
clear-cut on smaller spatial volumes, indicating the absence of (non-singular) finite volume corrections at A'^ = co 
- reminiscent of large A^ reduction. The observed rapid increase of the inter-phase surface tension with A' may 
indicate that for N = oo the deconfinement transition cannot, in practise, occur. 



1. INTRODUCTION 

Originally proposed as a possible way of solv- 
ing QCD (Ij, the large N limit of SU{N) gauge 
theories plays an important role for our under- 
standing of this class of theories. 
It has been shown numerically ^ that several 
quantities at zero temperature have a well-defined 
large N limit and that deviations from that limit 
at finite N are small and accounted for by a 
0{1/N'^) correction (the first expected correction 
by perturbative arguments) all the way down to 
N — 2. In [S] similar features were observed at 
finite temperature. In particular, for the decon- 
finement temperature Tc it was found 



Tc/V^ = 0.582(15) + 0.43(13)/A^2 



(1) 



Here we extend the investigation of [3 . We study 
SU{N) groups for iV = 2, 3, 4, 6, 8. New improved 
results for both Tc and ^/a' help us to keep under 
better control the extrapolation to infinite N. 
In the following we will discuss our numerical re- 
sults for the critical temperature and the latent 
heat Lh as a function of N and what the observed 
behaviour at finite A'^ may imply for the physics 
at finite temperature in the large N limit. This 
work is based on 0] , to which we refer the reader 
for a more detailed discussion. 



2. LATTICE SETUP 

In our numerical simulations we have used the 
Wilson action 



^ = /3E{l-^I^eTrC/p| 



(2) 
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where (3 — 2N/g'^, with g the coupling of the the- 
ory. Up is the path ordered product of the links 
around the plaquette P and the links are SU (N) 
matrices. The sum is performed over all the pla- 
quettes P. The lattice size is L in three directions 
and Nt (with Nt ^ L) in the fourth one. We call 
spatial directions the three of equal size and tem- 
poral direction the remaining one. 
The system can be regarded either as a statistical 
system in four dimensions with the role of tem- 
perature played by /3 or as a three-dimensional 
quantum field theory at finite temperature, with 
the physical temperature given by 

T, = l/Nta{P) , (3) 

a being the lattice spacing. Thcrmodynamical 
quantities have been obtained via a finite size 
study in the critical region of susceptibilities like 
the specific heat and the susceptibility of the 
Polyakov loop. 

3. CRITICAL TEMPERATURE 

The critical temperature is unambiguously de- 
fined only in the thermodynamic limit. On a fi- 
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nite lattice, one usually takes as a definition the 
value at which the susceptibility of some observ- 
able relevant for the transition has a peak. More 
specifically, one defines the pseudocritical cou- 
pling Pc{Nt, L) as the (3 at which the chosen ob- 
servable displays a peak at fixed L, Nt. Then the 
value of the critical coupling is given by a fit to 
the finite size scaling relation 

(3c{Nt,L = cx,)=p,{NuL) + hL-^/'' , (4) 

where, if d is the number of spatial dimensions, 
we get 1^=1/(1 for a first order phase transition 
and v < 1/d for a second order one. The value 
of l3c{Nt, L = oo) does not depend on the chosen 
observable. 

In principle the order of the phase transition can 
be extracted from the very same procedure used 
for determining Pc{L — oo). However, in practise 
that proves to be hard: using as a fit ansatz 
gives large errors on v. For this reason, we have 
determined this exponent from the scaling of the 
maximum of a susceptibility. In fact, from general 
arguments 



X7nax{Nt,L = oo) CX L'^^" 



(5) 



where 7 is the critical exponent governing the 
scaling of the considered observable. The fitted 
value must satisfy the constraint < j/v < d, 
with the upper bound saturated by first order 
phase transitions. 

The cleaner signal for the fit (0) is generally given 
by the specific heat, from which we get evidence 
for a first order phase transition for > 2. In- 
serting this information back into eq. we ob- 
tain /3c(-/Vt, 00). Next, by measuring the string 
tension at zero temperature at that value of (3, 
we get Tc/ \fa. 

We have determined /3c (i — 00) for at least three 

values of Nt^ chosen in such a way to avoid any 

bulk phase for TV > 4. This allows us to obtain 

the continuum limit of Tc/^/a'^. 

The final step is the extrapolation to = 00. We 

find 



Tc/V^ = 0.596(4) + 0.453(30)/7V2 



(6) 



^Our finite size analysis has been performed in full only at 
Nt = 5. For other values of Nt, we have studied one lattice 
size in the thermodynamic regime and then we have used 
the scaling with /3 of /i in eq. to determine /3c. 
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Figure 1. Extrapolation of Td \fo to infinite TV. 
The dashed line is the best fit to the data. 



Our data for Td -Ja as a function of l/N"^ are 
plotted in fig. ^ Despite the increased precision, 
which reflects on a noticeable reduction of the 
error bars with respect to still the leading 
correction in l/N"^ accurately describes the data 
all the way down to SU{2). 

4. LATENT HEAT 

At finite > 3 the deconfinement phase tran- 
sition is first order. First order phase transitions 
are characterised by a latent heat, i. e. a difference 
in action between the confined and the deconfined 
vacuum at T = Td. The latent heat (per site) 
can be obtained from the infinite volume limit of 
the maximum of the specific heat in the following 
way: 



LhiNt)^ lim 



PciL,Nt) 



C„tax{L,Nt) 



6i3 



(7) 



We have measured Lh for Nt — 5 and > 3, 
and we have found that its value increases with 
A^. Assuming a leading 0(1/A^^) correction, our 



^Here we are referring to the latent heat of the statistical 
system, which is related to the jump of the gluon con- 
densate across the transition in the quantum theoretical 
system. 
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Figure 2. Extrapolation of Lh to infinite N. The 
dashed curve is the best fit to the data. 



results at finite N can be extrapolated to iV = oo 
(cfr. fig. Ol- We obtain a finite value for Lh in 
that limit, which implies that the transition is 
first order at = oo. The negative coefficient 
for the correction indicates that the strength of 
the transition increases with N. 



N increases (a) the thermodynamical limit is 
reached earlier and (b) finite size corrections are 
smaller. This suggests that at large N the scaling 
regime is obtained for L = Nt + e, with e — > as 
N ^ oo. If we reduce L below Nt, the system 
undergoes a deconfinement phase transition due 
to the breaking of a spatial Z{N) symmetry. This 
has been investigated from a different perspective 
in The breaking of a spatial Z{N) symme- 
try could be prevented by twisting the boundary 
conditions in all spatial directions. In this way 
a finite temperature regime can be obtained for 
L < Nt. It is then plausible that at infinite N 
the spatial volume can be shrunk to a single pla- 
quette. This naturally connects to the twisted 
Eguchi-Kawai model (7j. 

6. CONCLUSIONS 

Our results for Tc and Lh indicates that a sensi- 
ble large N limit exists for SU{N) gauge theories 
also at finite temperature. Deviations from that 
limit at finite N are small and can be accounted 
for by a leading 0{1/N'^) correction. Our ex- 
tended range of N (up to eight) reveals some in- 
teresting scenarios for the physics at = cx3. At 
present we are extending this analysis to other 
observables. 



5. SURFACE AND FINITE SIZE EF- 
FECTS 

A consequence of the increasing strength of the 
phase transition is that a stronger interface sep- 
arates the confined and the deconfined vacua at 
larger N. Hence the life time of a metastable 
state at fixed spatial volume increases as N is 
increased, suggesting that it becomes infinite at 
TV = oo even for /3 close to but different from 
Pc- In that case, we would have two stable states 
in a finite region of /3's. Interpreted in terms of 
the Master Field 0, this suggests the possibility 
that this field is not unique: at infinite N there 
must be more than one Master Field, separated 
from each other by infinite energy barriers. The 
interplay between the different vacua at finite but 
large N guarantees a rich physics in the limiting 
case. 

When we vary the volume, we observe that as 
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